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ABSTRACT 

It has been well established that the /-mode of relativistic ordinary-fluid 
neutron stars displays a universal scaling behavior. Here we study whether the 
"ordinary" /q- and "superfluid" /g-modes of superfluid neutron stars also show 
similar universal behavior. We flrst consider a simple case where the neutron 
superfluid and normal fluid are decoupled, and with each fluid modeled by a 
polytropic equation of state. We flnd that the /o-mode obeys the same scaling 
laws as established for the /-mode of orindary-fluid stars. However, the oscillation 
frequency of the /g-mode obeys a different scaling law, which can be derived 
analytically from a homogenous two-fluid stellar model in Newtonian gravity. 
Next the coupling effect between the two fluids is studied via a parameterized 
model of entrainment. We flnd that the coupling in general breaks the universal 
behavior seen in the case of decoupled fluids. Based on a relativistic variational 
principle, an approximated expression is derived for the flrst-order shift of the 
/s-mode squared frequency due to the entrainment. 

Subject headings: dense matter — equation of state — gravitational waves — stars: 
neutron — stars: oscillations 



1. Introduction 

Ground-based gravitational wave detectors are either already operating or will soon 
be operating. For example, the LIGO detectors have conducted flve science runs since 
2002. While the current detectors are still not sensitive enough to detect gravitational waves 
directly, interesting upper limits on the gravitational- wave strains from several potential 
astrophysical sources (e.g., isolated pulsars and stochastic background) have been placed 
jAbbott et al.ll2007al lbll \t is very likely that the detection of gravitational waves will come 



within a decade or so after the upgrade of the detectors. 
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Pulsating neutron stars are one of the interesting potential sources of gravitational 
waves. They may also rotate rapidly enough for various kinds of ro tational induced mode 
instabilities to develop and enhance the emitted gravitational waves (lAnderssonI |2003| ) . The 
pulsation modes of neutron stars are damped due to the emission of gravitational waves, and 
hence are called quasinormal modes (see, e.g.. lKokkotas fc Schmidtlll999l : lFerrari fc Gualtieri 
20081 . for reviews), instead of normal modes as in Newtonian theory. The study of quasinor- 
mal modes of neutron stars has a long history dating back to the pioneering works of Thor ne 
and his collaborators (Thorne fc Campolattaro 1967 : Price fc Thorne 1969 : Thorne 1969al lbl: 



Campolattaro fc Thorne 197ol ). It is by now well established that the quasinormal mode 



spectrum of a neutron star is tremendously rich and the gravitational waves emitted from 
a pulsating star carry important information about th e internal structure of the star (e.g.. 



Andersson fc Kokkotas 1996, 1998: Benhar et al. 1999: Kokkotas et al. 2001: Benhar et al. 



2004 



Tsui fc Leung l2005al ). The gravitational- wave signals from pulsating neutron stars. 



should they be detected by ongoing or future detectors, can thus in principle be used to 
constrain the supranuclear equation of state (EOS), which is still poorly understood. 

While the quasinormal modes of a pulsating neutron star in general depend sensitively on 
the stellar mod el and EOS, some empirical un iversal behaviors have also been observed with 
different EOSs. I Andersson fc Kokkotas! (119981 ) noted that the frequencies and damping times 
of the leading gravitational- wave w-modes and fluid /-mode of nonrotating neutron stars can 
be approximated by some empiri cal relations which depend only on the mass and radius of 
the st ar for most EOSs (see also lBenhar et al.lll999l ). In particular, I Andersson &: Kokkotas 
( 119981 ) obtained the following formulas respectively for the real and imaginary parts of /- 
mode frequency: 

Re(cjM) = aiM + a2C^'^, (1) 

\m{ujM) = (AC + , (2) 

where ai, A, and /?2 are model-independent constants determined by curve fitting. M is 
the mass of the star and the dimensionless parameter C = M/R (in units where G = c = 1) 
is the compactness of the star. The real part of uj is the mode oscillation frequency, while 
the imaginary part is inversely proportional to the damping time due to the emission of 
gravitational waves. Recently, t he physical me c hanism behind such universal behaviors has 
been investigated in detail by iTsui fc Leung J2005bh . In their notations, the (complex) 
frequency uj of the if-mode or /-mode can be approximated by 



ujM = aC^ + bC + c, 



(3) 



where the complex parameters (a, b, c) are also model-indepe ndent constant s deter mined 
by curve fitting. Exploring such universality for the w-modes, iTsui fc LeungI (j2005al ) have 



- 3 - 



developed an inversion scheme to determine the mass, radius, density distribution, and the 
EOS of a neutron star model from the frequencies of the first few w-modes of the star. 

Mature neutron stars are known to be very cold on the nuclear temperature scale (10^'^ 
K) . It is believed that a newborn neutron star can cool down rather quickly (on a timescale 



Lombardo 


1999; 


Lombardo &; Schulzel 


2001 



Andersson et al.ll2005l ). For sufficiently high density, nuclear matter can transform to de- 
confi ned quark matter, which may also b e in the so-called color-superconducting phase (see. 



e.g 



Alford fc Reddyl l2003l : lAlfordll2004j ). To date, there is still no direct evidence for the 
existence of nucleon superfluidity in neutron stars. However, the well-established pulsar 
glitch phenomenon is best explain ed by the pinning and unpinning of large numb ers of su- 
perfluid vortices to the solid crust (IRadhakrishnan fc Manchesterlll969l : lLynelll993l ). So, how 
would nucleon superfluidity affect the dynamics and hence the oscillation mode spectrum 
of neutron stars? For the simplest model, neutron stars consist of neutrons, protons, and 
electrons. When neutrons become superfluid, they dynamically decouple from protons and 
electrons. The protons can also be in a superconducting state, but they are coupled to the 
ordinary electron fluid via electromagnetic intereaction on a very short timescale. In effect, 
the neutron star interior can be approximated by a two-fluid model: one fluid is the neutron 
superfluid; the other fluid is a conglomerate of all other charged constituents. 

To a flrst approximation, the two interpenetrating fluids could be considered as indepen- 
dent. However, due to the strong interaction between neutrons and protons, the two fluids in 
general can couple via the so-called entrainment effect. Entrainment is a multifluid effect and 
it arises when the flow of one fluid induces a momentum in the oth er fluid. In this work, w e 
will study the parameterized entrainment model used previously by I Andersson et al.l (120021 ). 



Building on the formalisrn developed by Carter and his collaborators (e.g . , 



Comer fc Langlois 



1993 



1994 



Carter fc Langloid Il995l ; iLanglois et al.l Il998l ) 



Cartel 


1989 


Comer et al. 



ing t wo-fluid stellar model (see, e.g., iLindblom fc Mendelllll994l ; lLeelll995l ; iPrix &: Rieutord 



2OO2I . for the corresponding study in Newtonian theory). These works show the existence 
of a new family of modes for a two-fluid star, the so-called superfluid modes. These modes 
have the distinguishing characteristic that the two fluids are essentially counter-moving, as 
a manifestation of an extra fluid degree of freedom which is missing in the study of pulsat- 
ing ordinary, single fluid neutron stars (jComer et al.lll999l ). Another characteristic of the 
superfl uid modes is that they depend sensitively on the entrainment effects between the two 
fluids (lAndersson et al.ll2002l ). In a two- fluid star, there is essentially a doubling of the fluid 
modes. For example, the single /-mode which exists in an ordinary fluid star is splitted into 
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an "ordinary" /-mode (denoted by /o), where the two fluids tend to move to gether, and a 



super fluid counterpart (denoted by /s) where the fluids are counter-moving (jComer et al. 



19991) ■ On the oth er hand, the ly-modes do not show such kind of mode doubling effect 



Comer et al.lll999l ). The two fluids of the w-modes always move in "lock-step" due to the 



fact that w-modes are largely spacetime oscillations. 

As mentioned above, the frequencies of the /- (or w-modes) of nonrotating ordinary-fluid 
neutron stars can be approximated by universal scaling laws ([I])- ([2]) for different EOSs. It is 
thus natural to ask whether this universality also holds for two-fluid neutron star models. In 
this paper, we investigate whether the /o and fs modes of two-fluid stars establish any kind 
of universality. In order to provide the reader a better understanding of the main results 
obtained in this paper, we shall give a brief summary in the following. 

In this work we first study the simplest case of two decoupled fluids, each with a poly- 
tropic EOS. We vary the polytropic indices of the two fluids to mimic the effects of different 
EOSs. We find that the /o-mode (i.e., the "ordinary" /-mode) can still be approximated 
very well by the same scaling laws ([I])-© as established for ordinary-fluid neutron stars. For 
the superfluid /s-mode, we find that the real part of the mode frequency satisfies a different 
universal scaling law (Eq. f|T4l) ). while the imaginary part in general does not have any uni- 
versal behavior. We are also able to derive the universal scaling curve for the real part of 
the /s-mode frequency analytically based on a homogeneous two-fluid model in Newtonian 
gravity. 

We then study the effect of coupling between the two fluids via a parameterized model 
of entrainment. We analyze both numerically and analytically how the entrainment affects 
the /s-mode frequency. Based on a relativistic variational principle, we have derived a 
general integral formula to calculate the first-order shift in the mode (squared) frequency 
(Eq. (1M|1 ). Furthermore, for the particular class of background EOS and entrainment models 
we considered in the study, we are able to approximate the integral formula by an algebraic 
expression (17^ which depends explicitly only on the model parameters. 

The main results obtained in this work (Eqs. f[T^ and flTTl) ) can be used by other 
researchers to obtain a good approximation to the /s-mode frequency (without the need of 
constructing their own numerical code) for the class of background EOS and entrainment 
models that have been used extensively to study superfluid neutron stars. 

The plan of this paper is as follows. Section [2] summarizes briefly the general relativistic 
two- fluid formalism. In section [3] we describe the EOS models we use in the study. Section H] 
presents our numerical results. Section [5] presents our analytical analysis for a homogeneous 
two-fluid star model. The effects of entrainment on the superfluid mode frequency is studied 
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analytically in Section [61 Finally, we summarize our results in section [71 We use units where 
G = c = 1 unless otherwise noted. 



The Two-Fluid Formalism 



In this paper we make use of the formalism and numerical code developed by I Comer et al. 



(I1999I ) to compute the oscillation modes of a general relativistic superfiuid star. As discussed 
in Sec. [H the superfiuid neutron star interior is approximated by a two-fluid model. One of 
the fluid is composed of superfiuid neutrons, and the other fluid contains all other charged 
particles (proton, electrons, and crust nuclei etc), which will be called proton fluid for simplic- 
ity. We focus our attention to a s implified mod e l whe re the two fluids are assumed t o coex ist 
throughout the whole star as in [Comer et al.l (Il999[ ). The work by [Comer et al.[ (Il999[ ) is 
based on Carter's gen eral relativistic superfiuid f ormalism. Here we will only summarize the 
formalism briefly (see [Andersson fc Comer[[2007l . for a recent review). 



The central quantity of the two-fluid formalism is the master function A(n^,p^ 



X 



which is formed by three scalars. 



n 



—na_n 



p 



—PaP°', and = —Uap'^- The four 



vectors n°' and are respectively the conserved number density currents for the neutrons 
and protons. The master function is a two-fluid analog of the EOS. Here we take —A to 
be the total thermodynamic energy density. Given a master function A, the stress-energy 
tensor is found to be 



(4) 



where is the generalized pressure 



<il = A-n''fx^-p''Xa- 



(5) 



The chemical potential covectors fia and Xa, respectively, for neutrons and (conglomerate) 
protons are given by 

fla = + Apa, Xa = + Au^, (6) 



where 



A 



OA 



-2^, C 



dA 



(7) 



It is noted that the so-called entrainment effect is described by the coefficient A. U A ^ 0, 
the mass current of one particle constituent will induce a momentum in the other constituent. 

The equations of motion for the two-fluid system consist of two conservation equations. 



V„n" = 0, VaP" 
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and two Euler equations, 



n"V[aHf3] = 0, p"V[«X/3] = 0. 



(9) 



Solving the linearized superfiuid and Einstein field equations, IComer et al.l (119991 ) obtained 
a system of first order differential equations to describe the (evei i-parity or po l ar) no n-radial 
oscillations of superfiuid neutron stars. We refer the readers to I Comer et al.l (119991 ) for the 
detailed equations and numerical techniques. 



3. Equation of State 



As discussed in the previous section, the master function A(n^,p^,x^) is a two- fluid 
analog of the EOS, replacing the role of P = P{p) (with P and p being respectively the 
pressure and total energy density) for ordinary one-fluid stellar models. While the univer- 
sality of the /-mode of ordinary fluid neutron stars is determined by many different EOSs, 
we cannot employ those EOSs directly for our two-fluid neutron stars. The reason is that 
realistic nuclear-matter EOSs are generally given in the form of pressure vs density. This 
is clearly insufficient for the two-fluid formalism used here in which the energy density —A 
must be given as a function of the two number densities (n, p), and henc e stand ard tabulated 
EOSs are incomplete for our purposes. We refer the reader to lLin et al.l (l2008l ) for a detailed 
discussion on what is required from next generation EOSs that allow for proper treatment 
of the dynamics of multi-fluid neutron stars. 

In this paper, we shall first study the oscillation modes of superfiuid neutron stars 
without considering the entrainment. The master function A is taken to be 

A = Ao = -m(n + p) - anu'^" - app'^", (10) 

where the two fluid particles are assumed to have the same mass m (i.e., the baryon mass). 



extensively to study superfiuid neutron stars ( 


Comer et al. 


I999I; 


Andersson & Comer 


2001b: 


Andersson et al. 


2002; 


Prix & Rieutord 


2002; 


Yoshida & Leel 


2OO3I). For this EOS, each fluid 



behaves as a relativistic polytrope and can be regarded as decoupled (though it should be 
noted that the fluids still couple globally through gravity). In general, we say that the 
two fluids are decoupled if the master function is separable, in the sense that it can be 
decomposed into two contributions corresponding respectively to each fluid, i.e., Kiji^.p^) = 
A„(n2)+Ap(p2). 
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3.1. Entrainment 

After studying the decoupled case, we shall also study the entrainment between the flu- 
ids. As already mentioned, entrainment describes essentially the effect where the momentum 
of one species depends on the relative motion between the two species. In the relativistic 
framework, this effect is specified by the coefficient A, which is in turn determined by the 
dependence of master function A on x^. Note that since the background i s assumed static. 



X — n p should remain small even when the star oscillates. We shall follow lAndersson et al. 



(120021 ) and approximate the master function by 

A = Ao + Ai(x^ -np), (11) 

where the first term Aq is given by Eq. (ITO!) . The second term is used to describe the 
entrainment and Ai is given by 

Al(n^p^) = e ^ y (12) 

T he positive constant e is used to parametrize the strength of the entrainment. As discussed 



m 



Andersson et al.l (120021 ). we shall consider the physically reasonable range to be < e < 
0.2. Note that the second term of the master function A vanishes in the background (where 
x"^ = np). Hence, the background quantities of the star are determined solely by Aq. We give 
the expressions for some of the thermodynamic coefficients on the background which will be 
used in later sections: 



A = -Xi{n\p'), 

B = --^-^A 

n on n 

C = -i^-^A (13) 
p op p 



4. Numerical Results 

4.1. Decoupled case 

In order to study whether the /o- and /g-modes of superfluid neutron stars establish 
any kind of universality, we have varied the values of (3^ and (3^ for the master function (Eq. 
(|T0|) ) to mimic the effects of different EOSs. The effects of and cTp on the mode frequencies 
are relatively small and hence we shall fix them to be o"n = o"p = 0.5m. Table [T] summarizes 
the different EOS models we shall study in this section. In Table [U we have defined two 
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parameters /3 = (/?n + /3p)/2 and A/5 = /3n — Pp. The parameter /? can be considered as a 
parameter that controls the bulk fluid motion, while A/5 represents the asymmetry between 
the two fluids. A larger value of A/5 would imply that the effects of the two-fluid dynamics 
become more important. On the other hand, the case A/5 = represents the one- fluid limit 
in which case the two fluid components are indistinguishable and the stellar models reduce 
to ordinary one-fluid neutron stars. 

In Figures [T] (a) and (b) we plot the real and imaginary parts of the /o-mode frequencies 
respectively against the compactness C = M/R for the EOS models Ai (with i = 1 — 5). 
These models have the same "bulk" EOS parameter (3 = 1.9, but different "asymmetry" 
parameters A/5. In particular. Model Al (with A/3 = 0) corresponds to the one-fluid limit 
for this series, while model A5 (with A/5 = 0.2) represents the most "asymmetric" one. For 
comparsion, we also plot the universal curves (Eqs. ([I])-([3])) for the /-mode of ordinary-fluid 
neutron stars in the figures. In Figured] (a), the solid line represents Eq. ([1]) for the star 
models with the masses and compactnesses obtained by EOS Al, while the dashed line 
represents Eq. ([3]). Similarly, in Figures [1] (b), the solid and dashed lines represent Eqs. ([2]) 
and ([3]) respectively. It is seen that the /o-modes of superfluid neutron stars still exhibit a 
universal behavior. For a given compactness C, the real and imaginary parts of the mode 
frequency do not depend sensitively on the EOS. Furthermore, the mode frequency can 
still be approximated very well by the universal curves ([I])- ([3]) satisfied by the /-mode of 
ordinary-fluid neutron stars. 

The corresponding results for the EOS models Bi and Ci (with i = 1 — 5) are plotted in 
Figures [2] and [31 It is seen clearly from Figures [1] to [3] that the /o-modes of superfluid neutron 
stars satisfy the same universal scaling as their ordinary- fluid counterpart. This result might 
not be so suprising as the fluid motions of the /o-mode are such that the two fluids move in 
"lock-step" and reach essentially the same maximum displacements at the stellar surface: a 
feature that is similar to the situation of an ordinary-fluid neutron star in which the (s ingle) 



fluid displacement reaches a maximum at the surface (see Fig. 9 of IComer et al.lll999l ). 



Now we consider the superfluid /s-mode, which does not have an ordinary-fluid counter- 
part. For this class of modes, there is no a prior leason to expect that they would establish 
a similar universal behavior. However, as shown in Figures HJ 0, and [HI the /g-modes do 
indeed follow another universal scaling law. In these figures, we plot the real parts of the 
/s-mode frequencies against the compactness C for the EOS models Ai, Bi, and Ci (with 
i = 2 — 5). In particular, the data can be approximated very well by the curve (dashed lines 
in the figures) 

Re{ujM) = {IC^Y^\ (14) 
where / is the spherical harmonics index. Note that the /g-mode does not exist in the models 
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Al, Bl, and CI since they correspond to the one- fluid hmit Aj3 = 0. 



In Figure [7] we focus on the model B5 and plot the real part of the /s-mode frequency 
against the compactness for different values of the spherical harmonics index / = 2, 3, 4. It is 
seen that the data can be approximated by Eq. f|T^ very well in general. The discrepancy 
between the numerical data and Eq. (fT4|) becomes significant only for the / = 4 modes in the 
region of high compactness {C > 0.2). In Sec. [5l we shall derive this scaling relation in the 
framework of Newtonian gravity and under the assumptions that the star is composed of two 
homogeneous and decoupled fluids. It should be noted that Eq. (HM does not contain any 
free parameters, while the scaling equations ([I])-© for the / -modes of ordinary-fluid stars 



depend on paramete rs that are determined by curve fitting (lAndersson fc Kokkotad Il998 
Tsui fc Leun"3l2005bh . 



It is also instructive to compare the scaling law ([T]) for the f-mo de of ordinary-fluid stars 
with E g. (fT4l) since they both contain the factor C^/^. As discussed by lAndersson fc Kokkotas 
(Il998l ). since the characteristic timescale is related to the mean density of the star p, it is 
expected that the /-mode frequency scales with p^/^ ~ (M /R'^Y^'^. Fitting to the n umer- 
ical data for the /-mode frequencies of ordinary-fluid stars, lAndersson fc Kokkotad (Il998l ) 
obtained the scaling law ([1]). While Eqs. ([1]) and f|T^ both contain the scale factor C^^^, we 
notice two important differences: (i) our scaling law f|T^ for the /g-modes does not contain 
any free parameter; (ii) it also does not contain a term that scales only with M as in Eq. ([1]). 
As we shall show in Sec. El the scaling law f[T^ is not a trivial generalization of Eq. ([T]) to 
the case of tw o-fluid stars. It is also not the well-known Kelvin mode for an incompressible 
stellar model (IChandrasekharl Il98ll ) . whose oscillation frequency is given by 



2l{l 



L (2/ + 1; 



1/2 



(15) 



However, Eqs. f[T4j) and f[T5|) agree with each other in the limit Z ^ 1. 



Now we turn to the imaginary part of the /g-mode. We see that in general it does not 
follow any universal scaling law. For a given compactness, its value depend very sensitively 
on the EOS models. In Figure [HI we plot Im(a;M) against C for models Bi {i = 2 — 5). 
These models have the same "bulk motion" parameter (3 = 2, but different "asymmetry" 
parameter A(3. It is seen that Im(u;M) decreases significantly with decreasing Af3. In fact, 
for the particular form of the master function equation (fTOl) . it is known th at the superfluid 
mode s become non-radiating, and hence lm{ujM) — * 0, in the limit A/5 — > (lAndersson et al. 
2002h . 



While Im(u;M) does not follow any universal law, it is interesting to note that Im(co'M) 
can be rescaled in such a way that the data from models with the same value of (3 lie on 
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the same curve. In Figure [9] we plot lm{ujM)/ Aj3'^ vs C for the models Ai, Bi, and Ci 
(i = 2 — 5). It is seen clearly that, for a given (3 (eg, the models Ai) and compactness C, the 
imaginary parts of the mode frequencies (normalized by the asymmetry parameter are 
rather insensitive to A(5. 

The above results suggest that, for the polytropic EOS model (fTOl) . the imaginary parts 
of the /g-modes are given approximately by 

Im(cjM) = /A/?2, (16) 

where the value of / depends only on jS and C. So why is there such a scaling law? In 
our models we have fixed o"n = cTp and thus we expect our equations to be invariant under 
A/? — > —A(5 and exchange of species label n and p. As lm{ujM) = when A/? = 0, we thus 
expect that the leading dependence of Im(ci;M) on A/5 must be quadratic. In fact , it ca n also 
be shown that the leading dependence is also quadratic even if 7^ CTp (IWongll2008l ). We 
shall present in the Appendix a qualitative argument to explain why the damping time of the 
/s-mode depends sensitively on various parameters in the EOS, including /3n, /3p, cxn and (Tp. 
Such dependence on the details of the EOS is definitely an advantage to asteroseismology 
as observation of /g-mode oscillations is then likely to provide clues to the coupling between 
neutron and proton fluids in neutron stars, which is not yet thoroughly understood. 



4.2. Effects of entrainment 

In the above we have demonstrated that the real part of the /g-mode exhibits a universal 
behavior. The value of Re(u;M) depends essentially only on the compactness C of the star 
(for a given value of /) and can be approximated very well by Eq. f[T^ . While we have 
only focused on a class of EOS given by Eq. ffTOj) . we believe that the universal behavior in 
general holds for decoupled fluids (as defined in Sec. [3]). However, it is also clear that one 
must consider the couplings between the two fluids in any realistic modeling of superfiuid 
neutron stars. Here we shall consider the entrainment between the two fluids using the 
master function ([TT]) in order to study the effects of coupling on the oscillation modes. It 
shou ld be noted that the e ntrainment effect on the ordinary /o-mode is rather small (see, 
e.g., Andersson et al.lbood ). Hence, we shall not consider the /o-mode in the study. 



We shall focus only on the real part of the /g-mode frequency and simply refer to it as 
the oscillation frequency hereafter. We use the two EOS models B2 and B5 (see Table [T]) 
for the first part Aq of the master function ffTTl) . It is recalled that the background stellar 
model is completely determined by Aq. The entrainment effect is given by the second part 
Ai and its strength is governed by the parameter e. For a given background EOS model Aq 
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and compactness C, we calculate the oscillation frequency of the /g-mode as a function of e. 
The results are summarized in Figures [TO] and [TTl 

In Figure[TD]we plot the squared /s-mode frequency i-o'(e)^ (normalized by its value when 
e = 0) against e for model B2. In the figure, the circle and square data points correspond to 
the compactness C = 0.1 and C = 0.2 respectively. The co rresponding results for model B5 



are shown in Figure [TTl As found by lAndersson et al.l (120021 ) . we see that the mode frequency 



increases with the entrainment parameter e. It is also seen that the rate of increase of the 
mode frequency depends on the compactness of the stellar models. 

It has been seen in Sec. 14.11 that the oscillation frequency of the /s-mode essentially 
does not depend on the chosen EOS model for a given compactness. Figures [10] and [TT] 
demonstrate that the entrainment between the two fluids in general breaks such kind of 
universal behavior. For the case where e = 0.1 and C = 0.1, the ratio a;(e)^/a;(0)^ is 1.193 
and 1.268 respectively for models B2 and B5. Note that the value uj{0) is the same for 
both models. Nevertheless, we shall show in Sec. 16.21 that the shift in the frequency for a 
given value of e, namely the ratio ct;(e)^/ct;(0)^, can be obtained anlytically from the following 
equation (see Eq. ([71])) which involves only the parameters of the background EOS model 
and the compactness: 



r- a \i/{;3n-i) 1 1 



coiOy"" ^ [1 + ^7(0.2)] 



(17) 



where the function g{x) is defined in Eq. ( I73|) . The excellent agreement between the analytic 
result and numerical data can be seen in Figures [10] and [TT] The solid and dashed lines in 
the figures represent Eq. (1171) respectively for C = 0.1 and 0.2. It is seen clearly that, for the 
physically reasonable range of e considered by us, Eq. (IT7|1 agrees with the numerical data 
very well. Significant derivation of the data from Eq. f[T7j) occurs only for higher values of e. 



5. Superfiuid mode for a homogeneous two- fluid star 

In this section, we shall derive the scaling law f[T^ for the real parts of the superfiuid 
/s-modes using the Newtonian two-fluid hydrodynamics equations. Our analysis is based 
on a model in which the densities of the two fluids are constant throughout the interior 
of the star. The general set of two-fluid equations in Newtonian gravity can be found in 
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Andersson fc Comerl (l2001al ): 








5tPp + 9i(ppfp), 

2a 

Pn 

2a 



Pn 



Pp 

2a 

+5*^9, ($ + pp) + —b'^buv'^d,{vl - vl) 
Pp 



Pp 



d'd,^ = 47rG(p„ + Pp), 



(18) 
(19) 

(20) 

(21) 
(22) 



where px, f*, and px (x = n, p) are the mass density, velocity, and chemical potential per unit 
mass of species x respectively. $ is the gravitational potential, a is a function to describe 
the entrainment and is defined by 



dU (pn, Pp, A^) = PnCipn + PpC^Pp + Otd^^ 



(23) 



where U is the internal energy density and = — t;pp. Note that the relativistic coeffi- 
cient A defined in Eq. ([7]) is proportional to a in the Newtonian limit (lAndersson fc Comer 



2nniah . 



For a non-rotating static background (f^ = 0), the equilibrium equations are given by 



<9^($ + Px) 



0, 



47rG(pn + Pp). 



(24) 
(25) 



In particular, we shall assume that the background star is homogeneous in which p^ and Pp 
are constants. The solutions for the background equations are then given by 



-¥(Pn + Pp)(3i?^-r^ 



-¥(Pn + Pp)i?^ 



f G(p„ + Pp) {R^ - r2) 



Px 



r < R 
r>R 

r < R 
r> R 



(26) 
(27) 
(28) 
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The linearized versions of Eqs . (fT8l)- (1221) have been stu died previously (e.g.. lLindblom &: Mendell 



I994J : lAndersson &: Comerll2001al : iPrix fc RieutordI |2002| ) . In particular, the linearized hy- 
drodynamics equations for a non-rotating static background are 



9f 



In 



Sp VSn Sp 



= -5,(5$ 
•^Pp + di (Pp4) = 0' 



5Pn), 



(29) 

(30) 

(31) 
(32) 
(33) 



where is the Lagrangian displacements for species x and S is used to denote Eulerian 
perturbations. For simplicity, we shall set a = in the following analysis. The effect of 
entrainment will be studied in Sec. O 

For our homogeneous background model, together with the assumption that the two 
fluids are decoupled, it can be shown from the linearized h ydrodynamics equations tha t the 
superfiuid modes are governed by the following equations (lAndersson fc Comeiil2001al ): 



+ V5(3 = 0, 



(34) 

(35) 



V • ^_ = 0, 

where d^_/dt = 6vp — 6Vj^ and 6P = 6fip — (5/in. A nonzero relative velocity ^_ is a 
characteristic of the superfiuid modes. For our simplified stellar model, we note that the 
equatio ns for the superfiuid rnodes are completely decoupled from those for the ordinary- fluid 
modes (lAndersson fc Comerl l2001al : iPrix fc RieutordI 120021 ) . Furthermore, the superfiuid 
modes are completely decoupled from the perturbation of the gravitational potential 5$. The 
counter-motion between the two fluids implies that the total density variation 6p = 6p^ + 6pp, 
and hence 5$, nearly vanish. Combining Eqs. (IMl) and (153]) . we obtain 

V^^/? = 0. (36) 

In solving Eq. (136|) . boundary conditions should be imposed at the center and the stellar 
surface: (i) the solution should be regular at r = 0; (ii) the Lagrangian variation oi (3 = 
Pp — Pn should vanish at r = i?. 



The background equation (12^ implies that pj^ = Pp, where a "prime'' 
derivative. The boundary condition at the surface can thus be written as 



denotes radial 



(37) 
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where is the radial component of the Lagrangian displacement of the species x. Eqs. (|3 
and fl37|) imply that 5(3 satisfies the following equation at the surface: 

^ - K^P' = 0. (38) 

Now let us decompose 5/5 in the form 5/3 = 6Pim{r)Yim{0, </))e*'^*. Eq. fl36|) and the regularity 
condition at r = imply that 

(5Am(r) = Ar\ (39) 



where A is some constant. Using the radial background profile for as given in Eq. (]28ll . 
the boundary condition (|38|) implies that the oscillation frequency uj of the superfluid modes 
is given by 

^ =^^- (40) 

This is the scaling law iHM presented in Sec. 14. H except that we have restored the gravi- 
tational constant G. As we have seen in Sec. 14.11 Eq. f[T^ can describe the data very well 
as long as C < 0.2. Relativistic effects become important for higher compactness. As can 
be seen from Figures [5] and [6l the derivation of the numerical data from Eq. pUj) becomes 
noticeable for C > 0.2. To further illustrate the relativistic corrections, we plot in Figure [12] 
the relative difference between the numerical data and the analytic result fHU]) for the fg- 
mode frequency of model C5. We see that the relative difference increases linearly with the 
compactness of the star. The relativistic correction rises up to about 8% for the maximum 
mass configuration. 

It is worthy to point out that Eq. fHU]) is in fact also the result for the /-mode of an 
incompressible single-fiuid star in the Cowling approximation (i.e., assuming 5$ = 0). For 
such a single-fiuid model, it can be shown that the perturbed one-fiuid equations governing 
the /-mode reduce to the same equation (136|) with Sp being replaced by the perturbed 
enthalpy Sh = 6P/ p: 

V^5h = 0. (41) 

The equivalence between the two governing equations fl5B]) and f HT]) comes from the fact that 
the counter-motion of a superfiuid mode leads to a nearly vanishing 5$ as discussed above. 
However, it should be noted that Eq. (HOj) does not describe the numerical data of ordinary 
/-modes as good as it does for the superfiuid /-modes. This suggests that incompressiblity 
and the Cowling appr oximation are not good approximations for the ordinary /-modes (see 



Andersson et al. 



20081 . for the analysis of the ordinary /-mode with the compressibility taken 



into account). 

Finally, if we also assume that the entrainment term a is a constant, Eq. fj40l) can be 
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easily generalized to include the effect of entrainment and becomes 



where the two parameters a and e are related by a = epnPp/[2(epn + (1 + e)Pp)]- At first 
sight, one might expect that this generalized result could be used directly to compare with 
the numerical data. However, it should be noticed that Eq. (H2l) involves the mass densities 
of the two fluid components. This makes the comparison not practical since the numerical 
data are obtained from stellar models where pn and pp are varying throughout the stars. 
In the next section, we shall use a variational principle to study the effect of entrainment 
perturbatively. 



6. Pertubative Analysis of the entrainment 

6.1. General integral formula 

In this section we shall study the effect of entrainment on the superfluid mode fre- 
quency for a two-fluid star without rotation. Our approach is to treat the entrainment as 
a perturbation and employ a variational principle to calculate the first-order shift in the 
mode frequency. We have generalized the v ariational principle for po lar oscillation modes of 
ordinary-fluid relativistic stars developed by lDetweiler &: Ipserl (Il973l ) to the case of two- fluid 
stars. The derivation of the general relativistic v ariational pr inciple is somewhat lengthy and 
we shall not present the details in this paper (see lWongl2008l ). However, in order to illustrate 
the basic idea, we shall first derive the corresponding variational principle in the Newtonian 
framework and then simply quote the relativistic result. 

The relevant equations for the Newtonian analysis are given by Eqs. (!29l) - (!33l) . We 
shall treat this system of equations as an eigenvalue system with the operator containing 
the entrainment term a as the perturbing Hamiltonian. Consider a normal-mode solution 
with the time dependence of the form (5/(r, t) = 6f{r)e^'^^, where 6f is a perturbed quantity. 
Eqs. f l29|) and fl30|) then become 



UJ 



UJ 



Pn 

. Pp 



(43) 



(44) 



We want to formally express the above equations (to first-order of smallness in a) in the 
form 

'^(0) + ^(1)) Y = u^Y, (45) 
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where is an operator i-th order in a and Y = ^p) is an abstract eigenvector with 
the two Lagrangian displacement vectors as its components. The system is solved with the 
"constraint" equations flMl) - fl5^ which do not contain time derivatives. 

First we define the "unperturbed Hamiltonian" ^(o) by 

^<o) ('') - ff'!*! " fr^] . (46) 



In order to apply standard perturbation theory as in quantum mechanics, we first need to 
define an inner product such that -ff(o) is symmetric. To this end, we define the inner product 

(^n, ^p) ■ (^n, ^p) = j dV{Pnri^ " In + Pp^p " Q- (47) 

We shall show that, under this inner product, H{q) is symmetric. In other words, {t}^^ rj^) ■ 
H(o) is symmetic in (t/^, r/p) and (ln,|p): 

(r7„,T7p)-^(o) 

(48) 

dSrj [pnrii{6^ + 6j2n) + Pp^^(5<I> + 6j2^ 



+ j dV SpniS^ + 5fl^) + Spp{S^ + Spp) 

where the first term is a surface integral which is carried out on the two sphere r = R and f 
is the outward unit normal to the sphere. To arrive at the second row, we have performed an 
integration by parts and used Eqs. fl3T|) and fl32|) . 6f and 6f refer to the Eulerian variation 
of the quantity / associated to the displacement and r}^ respectively. 

Note that the surface integral in Eq. (HHj) vanishes since the densities tend to zeros at the 
surface for our stellar models. The remaining volume integral can be shown to be symmetric 
in (r/jj,r7p) and (ln,|p). To show this, we first note that 



and thus 



d^U ^ dU , 

^Pn=^^Pn + ^5pp. (50) 
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It should be mentioned that a term associated with dA'^ vanishes because in the background 
A = 0. This imphes that 



dpi 

+ (SppSpa + Spjpp) 



p 
dU 



(51) 



dpndpp 



which is symmetric in (5pn, ^Pp) and ((5pn, ^Pp)- It can also be shown easily that / dV{6pn + 
(5pp)(5$ is symmetric. Thus the volume integral in Eg. (1481) . and hence -^(o), is symmetric. 

We can now apply standard perturbation theory as in quantum mechanics with the 
perturbing potential -ff(i) defined by 



^(1) 



2aiu 



(0) 



1/Pn -1/Pn 
-1/Pp 1/pp 



(52) 



where uj'^q^ denotes the squared frequency of the mode in the absence of entrainment (i.e. 
a = 0). The first-order shift in the squared frequency is given by 



UJ 



(1) 



(0) ^ (0)^ 
P ^ 



a 



(0) 



i (0) 
Sp 



(53) 



/^^(Pne^-e^+PpC^-lp' 



(0) . t (0) 



where and ^^'^^ are the Lagrangian displacement vectors associated to the "zeroth-order" 
mode solution with squared frequency u^^y 

For a given static background of two-fiuid stellar model, the effect of a small en- 
trainment on the oscillation mode frequency is to change the squared frequency u^q^ to 



cu^ = u^Q^+u'^^y It can also be seen qualitatively from Eq. (153|) that the ordinary-fluid modes 
in general do not depend on entrainment (as has been discussed in lAndersson fc Comer 
2001al : lAndersson et al.ll2002l : iPrix fc RieutordI l2002l ) . The first-order shift depends on the 
relative displacement |^n^°^ — ^p^°^|, which tends to zero for ordinary-fiuid modes. 

Eq. fl53l) can readily be used to compute the change in the mode frequency once the 
entrainment function a and the zeroth-order solution (^n^^^p^^) are given. However, we 
cannot employ the equation directly since it is derived in the Newtonian framework, while 
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our zeroth-order solution (^n"'', ^p'^'') are computed using a relativistic numerical code. Nev- 
ertheless, with the derivation of Eq. fl53l) as an illustration of the basic idea involved, we 
h ave derived the cor r espoii ding result based on a two-fluid formalism extension of the work 



of Detweiler fc 



lengthy ( Wongl 



pser (119731 ) ■ The derivation of the relativistic case is somewhat tedious and 



20081 ). Here we shall only present the final result of the first-order shift in 



LU 



4) 



(0) 



/ dVse-^/' Anp If - I 



(54) 



where dVs = e ' r sin OdrdOdcj) is the proper volume element for the spatial 3-geometry of 
the background spacetime metric 



ds^ = -e''^^)de + e^('^)rfr2 + [dO^ + sin^ edcj)^) . 



(55) 



The 3- vectors |f and ^p''-' are still defined as the Lagrangian displacement vectors associated 
to the zeroth-order solution as before. But the scalar product between any two 3-vectors U 
and V is now given by 



U V = e^VV + r'^U^V^ + sin^ OU'^V'^ , 



(56) 



where the ?7*'s and V^*'s are the coordinate components of the vectors. The entrainment is 
now described by the function A as discussed in Sec. [2J It can also be shown easily that 
Eq. fl54p reduces to Eq. fl53p in the Newtonian limit. 



6.2. Decoupled polytropes 

Eq. (15^ is a general result and could be used to obtain the frequency shift due to a 
small entrainment once the zeroth-order mode solution (i.e., without the entrainment) is 
given. In deriving Eq. (1541) we have not made any assumption about the master function 
A. In this subsection, we focus on the special case where the master function A is given by 
Eq. ffTTj) . We shall show that the integral formula (|5^ can be approximated very well by an 
algebraic relation which involves only the parameters of the background EOS (i.e., the part 
Aq in Eq. flTTl) ). the entrainment parameter e, and the compactness of the star C. 

To obtain an approximation to Eq. (IMI) . we first note that the chemical potentials of the 
neutrons and the (c onglomerate) protons evaluated on the static background are respectively 



fIComer et al.lll999f ) 



li = Bn + Ap = iioo^ X = Cp + A = XooG "^'^ , (57) 
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where fioo and Xoo are constants. The condition of chemical equihbrium, with Eq. f|T3l) for 
the thermodynamic coefficients, imphes that 



dn dp 

Using this relation, Eq. fl54l) can now be written as 



(5^ 



^(1) 



(0) 



jdV^e-"!^ {Anp \D\^) 



where the 3-vectors U and D are defined by 



(59) 



U 



n + p 



(60) 



We note that for the superfiuid modes where the two fiuids are dominated by counter- 
moving motion (i.e., when |-D| is large), Eq. (|59l) can be approximated by 



LU 



1 uj^^ J dV^e-^/^Anp \D\ 



(1) 



To show that this is in general a good approximation, we define the ratio 

^ / dV^e-'^ [{n + p) \U f] 



(61) 



(62) 



and present its value for the /g-modes of a few typical stellar models in Table El It is seen 
that 7] is in general much smaller than unity. Hence, it is a good approximation to neglect 
the first integral in the denominator of Eq. (l59|l . 



In the relativistic perturbative formalism (jComer et al.lll999l ). the Lagrangian displace- 
ments are decomposed into 

= e-^/V'-Wx(r)PKcos^), 

= -r'-V.(r)^PKcos^), (63) 
where Pi{x) is the Lengendre polynomial. Hence, Eq. (I6T!) can be put into the following form 



^1) 1 / drr'^'-e^'i" AnpD{i 



(64) 
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where 



+/(/ + l)[K(r)-V^p(r)]^ 
We also observe that the two functions (which are the two integrands in Eq. flM|) ^ 



(65) 



Mr) 
f2{r) 



r e 2 AnpDyr 
^21 A -u 



r 62 



n + p 



(66) 



have very similar shapes and both reach a maximum value at the position r ^ 0.8R (see 
Figure [T3|) . We thus further approximate Eq. fl6^ by 



UJ 



(1) 



(0) 



1 fijom 

/ioo /2(0.8/?) 

1 



n + p) ^ . (67) 

- r=0.8R 

For the decoupled polytropic model Aq (see Eq. ( ITOj) ) we use for the background, we have 



/iooC 2 = m + an/3nn' 



/XooC 2 = m + CTp/^pP^P ^ 

Since n(i?) = p{R) = and e^(^) = (1 - 2C), we deduce that 



Furthermore, it can be shown that 



mr) 



p{r) 



1 //Uoo 



Vl - 2C. 



(68) 
(69) 



(70) 



-e 2 



l/(/9p-l) 



(71) 
(72) 



where cXx = cTx/m, and 



m 



1 + 



c 



-X + 



C (2 - C) 



vl-2C)'^ ' 2(1-2C)' 



X 



(73) 



where x = 1 — r / R. 
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Using the above relations, we obtain the final expression 

(- a ^l/(/5n-l) 1_ 



^fo) (1 + ^7(0.2)) 



(74) 



It should be noted that this final formula only depends on the parameters of the master 
function and the compactness of the model. In particular, in contrast to the "exact" integral 
formula (Elj), it does not depend on the zeroth-order mode functions ^p^) which needed 
to be determined numerically. 

In Sec. 14.21 we have seen how well Eq. (1741) agrees to the numerical data. Here in 
Figure [14] we demonstrate the good agreement again by plotting the numerical data, the 
"exact" integral formula (15^ . and the approximated formula (JTlj) for the /g-modes of EOS 
model B2 at compactness C = 0.15. In the figure, we plot the (normalized) squared /s-mode 
frequency uj"^ against the entrainment parameter e (see Eq. (1121) ). The circle data points 
are obtained directly from the relativistic numerical code. The solid and dashed lines are 
obtained respectively from Eqs. ( l54l) and (l74l) . It is seen clearly the excellent agreement 
among the three results for the physical range of e we consider. 



7. Conclusions 

In this work we have studied whether the /o- and /s-modes of superfiuid neutron stars 
exhibit any kind of universal scaling laws as have been seen in the /-mode of ordinary-fiuid 
neutron star models. We first focus on the simplified case of two decoupled fiuids, each with 
a polytropic EOS. We vary the polytropic indices to mimic the effects of different EOSs. Our 
numerical results show that the /o-mode, where the two fluids move in "lock-step", obeys 
the same universal scaling laws as the /-mode of ordinary fluid stars. 

On the other hand, we flnd that the oscillation frequency of the /g-mode, which corre- 
sponds essentially to counter motion between the two fluids, obeys a different scaling law 
(Eq. f|T^ ). We have also derived the scaling law analytically based on a homogeneous two- 
fluid stellar model in Newtonian gravity. However, the damping time of the /g-mode in 
general does not exhibit any kind of universality. While we have only used a generalized 
polytropic EOS in our study, we believe that our conclusion holds in general for superfluid 
neutron star models in which the two fluids exist throughout the whole star and are de- 
coupled in the sense that the master function A (ie, the EOS) can be decomposed into two 
contributions corresponding to each fluid, i.e., A(n^,p^) = An(n^) -|- Ap(p^). 

The inclusion of a coupling term in the master function will in general break the universal 
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behavior. To illustrate the effect of coupling, we have studied the entrainment between the 
two fluids using a parameterized entrainment model. We show numerically how the /g- 
mode frequency increases with the strength of the entrainment. Furthermore, based on 
a relativistic variational principle, we have carried out a perturbative analysis and have 
derived an expression for the first-order shift of the frequency due to the entrainment. If 
the superfiuid /s-modes could be detected by future gravitational wave detectors, then the 
derivation of the observed mode frequencies from the universal scaling curve for the decoupled 
fluids could then be a useful probe to the coupling effects between the neutron superfiuid and 
normal fluids inside neutron stars. In summary, our main results (Eqs. f|T^ and (|T7|) ) can 
be used to obtain a good approximation to the oscillation frequency of the /g-modes for the 
generalized polytropic EO S and entrainment models that have been used ex t ensively to study 



superfiuid neutron stars (IComer et al. 



1999 



2OO2I : IPrix fc Rieutordlbood : IVoshida fc LeelboOsh T 



Andersson &: Comer 2001b: Andersson et al. 



While we have not yet detected the gravitational waves emitted by neutron stars, it is 
worthy to mention that the first gravitational-wave search sensitive to the /-modes has re- 
cently been carried out by the LIGO detectors and interesting upper limits on t he wave strain 



(with in the predicted range of some theoretical models) have also been placed (lAbbott et al. 



20081 ). In view of the fact that the advanced LIGO detectors will have more than a factor 
of 10 improvement on the sensitivity of the wave strain, gravitational-wave astroseismology 
may soon become a reality. 



This work is supported in part by the Hong Kong Research Grants Council (Grant No: 
401807) and the direct grant (Project ID: 2060330) from the Chinese University of Hong 
Kong. 



A. Gravitational-wave emission of superfiuid modes 



To gain some insight into the reason why the imaginary part of the /g-mode fails to 
follow a universal scaling which depends solely on the mass and radius of the star, we provide 
here a qulitative understanding based on the multipole f ormulas of the gravitational- wave 
luminosity for a single oscillation mode with frequency uj (jThornelll980l ): 



dE 
Itt 



(Al) 



1=2 



where Ni is some constant depending on the spherical harmonic index /. 6 Dim and 5 Jim are 
respectively the mass and current multipoles associated to the oscillation modes. For a non- 
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rotating two- fluid star with weak in ternal gravity (i.e., a Newtonian source), the multipoles 



are given by (lAndersson et al.ll2008l ) 



5Am = J{5p,, + 5p,yY;^dV, 

SJim = l^^Jr\pJv^ + p,6v^)-Yf^dV, (A2) 



where = [/(/ + 1)] f x Vyjm is the magnetic-type vector spherical harmonics. 

The ordinary /o-mode of a two-fluid star is characterized by the fact that the two fluids 
are comoving, which implies a large total density variation 6p = 6pn + Spp. Hence, similar 
to the /-mode of a single-fluid star, it is expected that the gravitational-wave emission of a 
/o-mode is dominated by the mass multipole SDim- The damping time of the / = 2 /-mode 
can be estimated by r ~ Em/uj^\SD22\'^- The mode energy Em can be calculated by giving 
the eigenfunction of the mode. For a homogenous single-fluid model, it c an be shown tha t 



the damping timescale depends only on the mass and radius of the star (IDetweileii Il975l ). 
In fact, the proposal for the leading scaling term in Eq. ([2]) was based on this rough 
estimation. Since the /o-mode of a two-fluid star is essentially the same as the standard 
/-mode, it is thus not surprising that the /o-mode also follows the same universal scaling 
law as we have seen in the numerical data. 

Now let us turn to the superfluid /g-mode. This mode is characterized by the counter- 
moving motion of the two fluids in such a way that the total density variation, and hence 
the mass multipole, nearly vanish. It is thus conceivable that the current multipole 6 Jim 
could provide the main radiation mechanism for the /g-mode. It is also interesting to note 
that, for the case of two nearly symmetric fluids (e.g., ~ dp and P^i ~ Pp in the master 
function (ITUl) ). the superfluid modes become non-radiating as the current multipole nearly 
vanishes. This explains why the imaginary part of the mode complex frequency Im(co'M) 
in the limit Af3 ^ as shown in Figure ([8]). In the general situation, however, the mass 
current PnSvn + Pp^Vp would depend on the mass fractions of the two components, which 
are determined by the condition of chemical equilibrium. The equilibrium condition in turn 
depends sensitively on the underlying EOS models and thus the damping timescale of the 
/s-mode would not have a simple scaling relation with the compactness of the star. Stellar 
models with the same global parameters (e.g., the compactness) could have vastly different 
fractions of the two particle species. The counter-moving character of the superfluid /s-mode 
would then lead to different amount of gravitational- wave emission among the stellar models. 

Finally, we note that it might be possible to perform a more detailed quantitative 
analysis by expanding the current multipole and extract the leading dependence of the 
damping time of the /s-mode on the local thermodynamics quantities. Such kind of analysis 
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might explain the relation f|T6|) seen in t he numerical dat a . A n atural starting point would 
be the extension of the recent work of lAndersson et al.l (120081 ) . in which the dependence 
of the damping time of the ordinary /o-mode due to the so-called mutual friction on the 
thermodynamics quantities has been studied. 
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Table 1. Models for the "polytropic" EOS defined in equation ffTOj) 



ModeP 




Pp 


oh 
P 




Al 


1.9 


1.9 


1.9 


0.0 


A2 


1.905 


1.895 


1.9 


0.01 


A3 


1.925 


1.875 


1.9 


0.05 


A4 


1.95 


1.85 


1.9 


0.1 


A5 


2.0 


1.8 


1.9 


0.2 


Bl 


2.0 


2.0 


2.0 


0.0 


B2 


2.005 


1.995 


2.0 


0.01 


B3 


2.025 


1.975 


2.0 


0.05 


B4 


2.05 


1.95 


2.0 


0.1 


B5 


2.1 


1.9 


2.0 


0.2 


CI 


2.1 


2.1 


2.1 


0.0 


C2 


2.105 


2.095 


2.1 


0.01 


C3 


2.125 


2.075 


2.1 


0.05 


C4 


2.15 


2.05 


2.1 


0.1 


C5 


2.2 


2.0 


2.1 


0.2 


^The paramet 


ers (Tn = 


(jp = 


0.5m are fixed in all models. 



b^=(/3„ + /?p)/2. 
'Af3 = (3n- /3p. 



Table 2: Values of the ratio rj as defined in Eq. fl62|) for the /g-modes of models B2 and B5 
(see Tabled]) at various compactness C 



Model 


C 


V 




0.10 


1.89 X 10" 


-5 


B2 


0.15 


2.24 X 10- 


-5 




0.20 


2.83 X 10' 


-5 




0.10 


6.95 X 10- 


-3 


B5 


0.15 


8.39 X 10- 


-3 




0.20 


1.11 X 10- 


-2 
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Fig. 1. — (a) The real and (b) the imaginary parts of cuM of the /o-mode are plotted against 
the compactness C for neutron stars described by the EOS models Az (z = 1 — 5). The solid 
line in (a) represents Eq. ([1]) for stars models with M and C obtained by EOS Al, while the 
solid line in (b) represents Eq. ([2]). The dashed lines in both figures represent Eq. ([3]). 
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Fig. 2. — Similar to Figured] but for the EOS models Bi (i = 1 — 5). The solid line in (a) 
represents Eq. ([1]) for stars models with M and C obtained by EOS Bl, while the solid line 
in (b) represents Eq. ([2]). The dashed lines in both figures represent Eq. ([3]). 
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Fig. 3. — Similar to Figure [H but for the EOS models Ci (z = 1 — 5). The solid line in (a) 
represents Eq. ([1]) for stars models with M and C obtained by EOS CI, while the solid line 
in (b) represents Eq. ([2]). The dashed lines in both figures represent Eq. ([3]). 
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Fig. 4. — The real part of uoM of the /g-mode are plotted against the compactnesss C for 
neutron stars described by the EOS models hi (z = 2 — 5). The dashed line represents the 
analytic result (Eq. (fT^ ). 
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Fig. 5. — Similar to Figure H] but for the models Bi (i = 2 — 5). 
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Fig. 6. — Similar to Figure H] but for the models Ci (i = 2 — 5). 
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Fig. 7. — Effects of the spherical harmonics index /: Re(to'M) of the /g-modes vs C for 
neutron stars described by the EOS model B5. The lines represent Eq. for different 
values of /. 
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Fig. 8. — Im(u;M) of the /g-modes vs C for the EOS models Bi (with z = 2 — 5). 
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Fig. 9. — Im(c<jM)/A/3^ of the /g-modes vs C for the EOS models Ai, Bi, and Gi (with 
i = 2-5). 
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Fig. 10. — /s-mode squared frequency co'(e)^ normalised by uj{Oy (the squared frequency 
when there is no entrainment) versus e for model B2 at various compactness C. The circle 
and square symbols represent the numerical data. The solid and dashed lines represent the 
results of Eq. (fTTD. 
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Fig. 11. — Similar to Figure fTOl but for model B5. 
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Fig. 12. — Relative difference (IcUnum — t^anal/t^ana) between the /s-mode frequency obtained 
by the relativistic numerical code tUnum and the Newtonian analytic result fHUl) c^ana versus 
the compactness C for model C5. 
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E 

Fig. 14. — /s-mode squared frequency uj{ey normalised by (^{0)"^ (i.e., its value when there 
is no entrainment) versus e for model B2 at compactness C = 0.15. The circle data points 
are obtained directly from the relativistic numerical code. The solid and dashed lines are 
obtained respectively from Eqs. ( l54l) and f l74l) . 



